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Abstract 

We investigate the magnetic property of the cold quark matter employing the 
magnetic moment p! as the order parameter. Through analysis of the effective po- 
tential V(p 2 ), we find that, at relatively high densities p > p c , the quark matter is in 
the normal phase. At p — p c , the magnetic phase transition takes place and, in the 
low-density region p < p c , the quark matter is in the ferromagnetic phase. For up, 
down and strange quarks, p c = 1.62p , 1.47p and 2.20p (Po the nuclear density), 
respectively. We also find that the leptons (e and p) inside the quark matter are in 
the normal phase. 
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1 Introduction 



Several years ago, the possibility of existence of a quark matter in a ferromagnetic 
phase has been pointed out [1]. Meanwhile, the possible discovery of a quark star 
[2, 3] has renewed our interest in this issue. 

In Ref. 1), through computing the energy density of a polarized quark matter, it is 
suggested that the Hartree-Fock state shows a spontaneous magnetic instability at low 
densities. The QCD coupling constant taken in Ref. 1) is a s (= g 2 / (47r)) = 2.2 and the 
quark mass is m = 300 MeV/c 2 . In this paper, we examine a possibility of the quark 
matter in a ferromagnetic phase by choosing a magnetic moment fl' oc (e/mtyirip as 
an order parameter. 

The QCD phase diagram is sketched in Fig. 1, where 'T' is the temperature and 
'//b' is the baryonic chemical potential. 'QGP' stands for the quark-gluon plasma 
phase, 'Hadron' stands for the hadronic phase and 'CSC stands for several phases 
which include the color super conducting phases. We are interested in the QGP phase 
near the phase boundaries and 0.6 GeV < In this region T << so that we 
ignore the effect due to finite T. Thus, the object of our concern is the quark matter 
with finite baryon-number density at zero temperature. 

We have to introduce two kinds of cutoff and then the result depends on the ratio 
of these two. This is not particular to the present study but is a common feature in 
the study of effective potential for some order parameter (see, e.g., Ref. 4)). We shall 
simply equate the two cutoff parameters (Sec. 3), for which we take the baryonic 
chemical potential \ib (Sec. 5). 

The plan of the paper is as follows. We introduce in Sec. 2 the closed-time-path 
(CTP) formalism [5, 6, 7, 8] and derive the formula for the expectation value (jl') 

— * 

under the existence of an uniform external magnetic field B. In Sec. 3, within the 
leading logarithmic approximation, we compute (fl') up to 0{B 2 B) for B and to 0(g 2 ) 
for the QCD coupling constant. Then, performing the resummation of the leading 
contributions, we obtain the "resummed form" for (fl'). In Sec. 4, we deal with the 
effective potential V((fl')). In Sec. 5, numerical analysis of V((fl')) is made in the 
range 0.6 GeV < fi B < 4 GeV. We find that, in the high-density region fi B > /i Bc , 
the quark matter is in a normal phase. As the density is lowered, the quark matter 
undergoes a magnetic phase transition, at the critical density, into a ferromagnetic 
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phase. The critical or transition chemical potentials are hb c = 641, 625 and 693 MeV 
for up (u), down (d) and strange (s) quarks, respectively. We finally discuss the case 
of leptons and find that they are in the normal phase. Sec. 6 is devoted to summary 
and outlook. 

2 Preliminary 

Throughout this paper, we adopt the leading logarithmic approximation, 

K<ln(A/m), 1« m[// B /(3m)] . (1) 

Here, A is the cutoff parameter and hb is the chemical potential being conjugate to 
the baryon number. The terms whose relative order of magnitude are 0(1/ ln(A/m)) 
and/or 0(1/ ln[/i B /(3m)]) are ignored. 

As the order parameter, we adopt one-half of the spin magnetic moment, 

$ = ^ = \ (^) ^ = ^757o# , 

where and in the following the color index is suppressed. An argument for this choice 
is given in Appendix A. 

2.1 Closed-time-path formalism 

We use the closed-time-path (CTP) formalism [5, 6, 7, 8], which is formulated by 
introducing an oriented closed-time path 7\ <g> T 2 in a complex time plane, that goes 
from — oo to +oo (71) and then returned from +oo to — oo (T 2 ). The field <p{x) with 
x eTi (i = 1, 2) is called type i, 4>i(x) = <p(x) I . The Lagrangian density in this 
formalism is t (= £(0i) — £(0 2 ))- 

The Lagrangian of a quark with some flavor reads 
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where £qce> is the interaction part of £qcd- It is convenient to construct a (2 x 2)- 
matrix "propagator" G, the inverse of the kernel of Eq. (2), which is diagonal in a 
color space. For the purpose of later use (cf. Eq. (8)), it is sufficient to obtain G with 
Bi = B 2 (= B). Straightforward manipulation yields for the (i,j) -element of G, 



G n (P) = -G* 22 {P)=G F {P)-N{p Q )[G F {P)-G F {P)\, 
G 12 (P) = [6(-p )-N(p )][G F (P)-G F (P)], 
G 21 (P) = [e( Po )-N( Po )}[G F (P)-G* F (P)}, 



(3) 
(4) 
(5) 



where 



Gf(P) 



A F (P) 



+ m + ^757o7 • B 



Am^ 



■— 2-7 • B + m7 5 7o - 75/70 
Am 



+2^ Po (-^B 2 l0 + m 75 7 • B - 75^7 • B 
Am \ Am 



x 



A ^ (P) - 2 (5^) 2(A/(P) + 2 ^ 2 



(6) 



P 2 -m 2 + i0+ 

Here we have used a capital letter for denoting a four vector P = (po,p) (p = \p\) and 
E p = ^/p 1 + m 2 . In the above equations, " * " does not apply to the Dirac gamma 
matrices and N(p°) is the number density function, 



N(po) = O(p o )0( f jiB/Z-Po). 



(7) 



Generating functional and the formula for (ft') 



The CTP generating functional is defined [7] by 



2 

Z(B U B 2 ) = J ]J [v^V^VAj] exp % J d 4 x (£ QCD + A ' #i - & ■ B 2 ) 

j'=i 



where Aj (j = 1,2) is the gluon field and p is the density matrix. An expectation 
value of the magnetic moment fl' is computed through [7] 



B 1 =B 2 =B 



I 

2 



JB! 5B 2 



\nZ(B 1 ,B 2 ) 



Bi=B 2 =B 



(8) 



We compute (fl') up to 0(B 2 B). 



3 Computation of (jj,') 



3.1 Oth order contribution 

From Eq. (8), we obtain for the Oth-order contribution (£qcd = 0) to (fl'), 



(?) = 



i e q _ Ti f d 4 P 



2 Am 



(2tt)< 



[7s7o7 (G n (P) + G 22 (P))} 



= ^ 7 ^/S [757ofGii(p)] - 

Straightforward manipulation using Eq. (6) yields, up to 0(B 2 B), 

d A P 



ImTr 



Ae q B T f d 4 P 
-Im 



3 Am"" J (2tt) 4 
4 / e \ 3 



(^1.) B 2 BIm 
15 \AmJ 



3A F + 2(2p 2 + 3m 2 )A 2 F 
d 4 P 



(2k) 



15A 2 



+A0(2p 2 + 3m 2 )A F + 8(8p 4 + 20m 2 p 2 + 15m 4 )A^ 



(9) 



(10) 



Substituting Eq. (3) for G n in Eq. (9) and using Eq. (10), one can compute (fx'). 
We divide (fx') into two pieces, 

(fl') = (/2')vac + {0)mat , 

where (fx') vac ((/?') mat) stands for the contribution from the vacuum (matter) sector. 
(jx') yac is given by Eq. (9) with G F for G u (cf Eq. (3)). Thanks to the 0(3,1) 
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symmetry, in Eq. (10), the replacements, p 2 — > — 3P 2 /4 and p 4 — > — 5(P 2 ) 2 /8 may 
be made. Thus we obtain 

in) 

In the second term on the right-hand side (RHS), we have dropped the terms that are 
proportional to m 2 (e ? /4m) 4 and m 4 (e g /4m) 4 , the contribution of which is nonleading 
in our approximation (cf. Eq. (1)). Introducing a momentum cutoff A, we have 

/S AMP) = ^ ln(A/m) ' (12) 

with which we obtain the first terms on the RHS's of Eqs. (28) and (29) below. 
For {Jl') mat , we obtain 



x 



(27T) 

3A F + VA| + 1 ( e -f- ) (lSA 2 , + 80p 2 A 3 F + 6 4p 4 A 4 F ) 



5 \ 4m 

(13) 

Here the terms that are proportional to m 2 (e q /Am) 2n and m 4 (e,/4m) 2 " (n = 1,2) 
have been dropped, which yield the subleading contributions. 

3.2 0(g 2 ) QCD correction to (//') 

Here we compute the 0(g 2 ) QCD correction to (//') in Eq. (8). Gluon propagator is 
diagonal in a color space. Choosing a covariant gauge, we have 

mo) = -(mQ)r = mQ) = -^ + vj^y a , w 

D^(Q) = 0(-q O )[D7(Q)-(D»/(Q))*}, (15) 

AT(Q) = WW]- (16) 

To 0(g 2 ), we obtain from Eq. (8), 

W = ? Tr / (S 75707 (K^^^)} + {i - 2 } ] , 



Ei,(P) = j9 2 (-) j+l J ^rGAP + Qh'D^Q). 
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Using Eqs. (3) - (5) and (14) - (16), we get, after some manipulation, 

W = ^rnTrf ^757o7[l - 2N{p )]G F {P)t F (P)G F (P) , (17) 

t F (p) = t n (P) + e( Po )t 12 (p) + e(-p )t 21 (p) 

= j9 2 J ^ 4 l fl [Gn(P + Q)D7(Q)-9( Po )G 12 (P + Q)D^(Q) 

-e(-p )G 21 (P + Q)D^(Q)]Y 
= jg 2 f^^lG n (P + QWAQ) 

+N( Po + g ) {G F (P + Q)- G* F (P + Q)} 

x {e(po)D%(Q) + e(-po)D^(Q)}} Y ■ 

The term with 0(—p°)D\ 2 (Q) in the last line vanishes, since D\ 2 (Q) = for g > 0, 
N(p + q ) = forpo + <?o < 0. Nevertheless we include it for later convenience. Using 
Eqs. (3), (15), and (16), we obtain 

Ai 2 r d^Q 



z F (P) = {g 2 Jj^YiG F (P + QWAQ) 



(2k)< 

-N( Po + g ) {G F (P + Q)- G* F (P + Q)} 
x {e( Po ) (0(-q o )D v /(Q) + 0(q ) W(Q))*) 
+0(- Po ) (6(q )D»/(Q) + 9(-q ) W(Q))*)}j T 

Substituting this for f] F (P) in Eq. (17), we finally obtain 

ae) = {/r)i + (/r) 2 , (is) 



( ^ = 3 9 A^ ReTr J (2^7 (2^ 757 ° 7 



x [1 - 2iV(p )]G' F (P)7' 1 G' F (P + Q)D F a {Q) 1 v G F {P) , (19) 
(A* /2 = -77^ T~ ReTr / 77m / 77m757o7 



r Am J (27r) 4 7 (2vr) z 
x [1 - 2iV(p )]G' F (P)7' 1 iV(p + g ) {G F (P + Q) - G* F {P + Q)} 

x {e(poW/(Q) + ^(-po)P>^(g)}7^ F (P) , (20) 
where D v ^ (D^ 1 ) is the retarded (advanced) gluon propagator; 

D R (Q) = [D A (Q)}* = 9(q )D F (Q) + 9(-q )[D F (Q)}* . 
The vacuum-sector contribution is given by Eq. (19) with N(p ) = 0. 
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It is necessary to compute (up to 0(B 2 B)) the quantity, 

Q = Tr h^oiG F (P)rG F (P'h,G F (P)] D T (Q) (P' = P + Q) , 
where T stands for F, R, or A. Introducing 

iZl = [A F (P)] l [A F (PT l ([A T (Q)] m=o r (T = F, R, A) , 
we obtain using Eq. (6), 



c (i) 

»self 



?(1)' 



C (2) 
»self 



4m 

-8 



4Vi + {3(g 2 + 2m 2 ) - tp. q + 4p g } 42 



+V(Q 2 + 2m 2 )xS 



3 

32 2 



3p p' -p-p' 



r {T) 
-221 i 



-16 



11X (T) 



-HQ- + bp 
4 



'411 



-8 [-9 + Xjg - JS) + 8 (2p Q p' - p 2 -p 2 -\p- p') X^ 



(21) 

(22) 
(23) 

(24) 



(T) 
221 



+8p' 2 X (T) 



230 



32 
15 



SpoPoP — 3p p —p p-p' — (p-p') )X 



r (T) 
'231 



(T) 



x 320 



-8p' 2 xS - 8 P 2 XS + ^ (3pV 2 + (P • P'f) X 33 o 
+ yP 2j 32o + ^ (-3pV 2 - p 2 p- p' + 5p p' p 2 - (p • p') 2 ) J j 



(T) 
321 



(25) 

Here C/^J ^ comes from the diagram that includes a quark self-energy subdiagram. 
For C/ggi f + £( 2 )', the contributions that are proportional to m 2n (n > 1) has been 
dropped, which do not yield leading contributions. 

According to Eqs. (21) - (25), we write (j = 1,2) in Eq. (18), with obvious 
notation, as 

= (£'>f If) + Wj (J = l,2). (26) 
The quark self-energy part £(P) consists of two pieces, £(P) = E vac (P) + E mat (P). 
Here E vac (P) stands for the self-energy part in vacuum theory, while E mat (P) is due 
to the existence of matter. Eqs. (19) and (20) tell us that (fl')^ 1 ^ includes E vac (P) 
and (/Z') 2 se1 ^ includes E mat (P). 



Observation on Eqs. (19) and (20) 

It is not difficult to see that, both in Eq. (19) and in (20), an integration over P 
(with Q fixed) is free from ultraviolet divergence. Then, we introduce a Pauli-Villars 
regulator (with mass Apy) for a gluon propagator, which makes the integrals in Eq. 
(19) and in Eq. (20) converge. This means that the contribution from the gauge- 
parameter (77) dependent part of the gluon propagator (Eq. (14)) vanishes, thanks to 
the Ward identity, and the results (28) - (36) below are independent of 77. Thus, we 
use the Feynman gauge, rj = 0. In the following we equate A PV with A in Eq. (12). 

3.3 Form for (//') 

Here we display the resultant form for (jl f ) , 



(//>= Ci + B 2 ^ 



B. 



(27) 



Summing over colors, we have for the contributions to C\ and to Ci from the vacuum 
sector and those owing to the existence of matter, in respective order, 



and 



(vac) 



£i(vac) 



(mat) 



£<(mat) 



3e 



32tt 2 

„4 



1 , 4a s A 
1 + — In — 

Sir m 



In-, 

m 



5127T 2 



rrr 



1 - 



4a, . A 
In — 

3ir m 



In 



A 



m 



H /4 

32tt 2 9m 2 
^4 -, 



512tt 2 m 4 



Aa s , 3A 
m — 



3tt 



4a, , 3A 2 
In 



3n m/iB 



ln^ 

3m 



(28) 
(29) 

(30) 
(31) 



Note that C\ 



(vac) 



(vac) 



and C[ 



(mat) 



are "exact" (cf. Eqs. (11) and 
(12)). To logarithmic accuracy, ln/xg/(3m) ~ ln/x^/m. Nevertheless we choose 
In li B /(3m) because the baryon number of the quark is 1/3 (cf. Eq. (7)). Adding 
both contributions, we obtain 

4a, 



c{ 0) = 



c 



(vac) 



c 



(mat) 



3e 



■ In — 

32tt 2 V m 



cf + 
/4 \ 



r* 1 



(i) 



(J = l,2), 



9m 2 



(32) 
(33) 
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) 



9m 2 hb 



2 



/4 3A 



(34) 



(35) 



(36) 



Several comments on the 0(g 2 ) contributions are in order. To logarithmic accu- 
racy, the contributions (/2') 2 (Eq. (20)) and (£[')[ (Eq. (26)) are negligible. In other 
ward, the leading contributions displayed above come from (jl')f eli \ to which the 
responsible diagrams include the vacuum piece S vac of the quark self-energy part. We 
see the point in more detail. 

On c| vac ) g (pi')'i- As mentioned above, integration over P (with Q fixed) in Eq. 
(19) converges. If the P-integration yields the term that is proportional to 
ln(\/— Q' 2 /m), the remaining Q-integration yields ln 2 (A/m) term. This is, how- 
ever, not the case. 

On G (pf)i' P-integration (with Q fixed) in Eq. (19) yields the terms being 

proportional to ln(v / — Q 2 /m), which lead to the ln 2 (A/m) terms. Cancellation 
occurs among them and no ln 2 (A/m) term remains. 

On C , 2 inat ' ) G (fl')'i'- Q- integration using the regularized gluon propagator yields 
ln(A/ra) terms. Cancellation occurs among those terms and then no ln(A/ra) 
term remains. If the Q-integration yields the ln[(m 2 — P 2 )/m 2 ] term, the re- 
maining P-integration results in the leading ln 2 (/is/(3m)) contribution. Among 
the Q-integrated terms are the ln[(m 2 — P 2 )/m 2 ] ones. Adding all of those, can- 
cellation occurs again, so that no ln 2 (yU B /(3m)) term remains. 

On {[!') 2 : Cancellation occurs between the leading contributions, so that there re- 



3.4 Resummation 

As has been mentioned above, the leading contributions (28) - (31) come from (/2')i Self \ 
for which the diagrams include S vac . In the region (4a s /37r) In A/m ~ 1 and/or 



mains no leading contribution in C[ vaiC \ C 1 mat \ and C^^. 
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(4a s /37r) ln// B /(3m) ~ 1, the expressions (28) - (31) are not reliable and an improve- 
ment is necessary through resummation of S vac - To this end, it is convenient to start 
from Eqs. (11) and (13) and, on the RHS's of them, make the substitutions 



P -> P(l - = Z?P(1 - Y?J , 
m -> m + S m = Z 2 -1 Z m (m + E^J . 



(37) 



Here E^ = Tr(/E vac )/(4P 2 ), E m = TrE vac /4 and Z 2 (Z m ) is the wave-function (mass) 
renormalization constant, 

7 2a A 

Z 2 = 1 - — In — , 

Sir m 

2a s , A 
Z m = 1 + — In-. 

7r m 

The substitutions (37) are not made for m within e q /(4m) and for p in 9(/i B /3 — po) 
in Eq. (13). It should be noted that we are using the Feynman gauge here. 

The leading-logarithmic contributions (28) - (31) emerge from Z 2 , Z rn and the 
piece ln(m 2 — P 2 )/m 2 in E^ and T! m . Then, to logarithmic accuracy, it is sufficient 
to resum the ln(m 2 — P 2 )/m 2 pieces in T/ w and Y! m , which is gauge-parameter (77) 
independent. Thus, Eq. (37) may be simplified as 



Z 2 l P 1 



3tt 



In 



P 2 



m 1 - 



m 



Zn 1 Z m m 1 



Aa s , m 2 - P 2 
In ^ — 

3tt m 2 



(38) 
(39) 



Straightforward manipulation yields 



(vac) 



^-,(vac) 



(mat) 



£,(mat) 



3e 2 l-5^1nA + 7(^l n A 

q y2 y2 37r m ' \ 3w rn 



fl_^l n AV 



1 A 

In — , 

m 



1024vr m 4 a. 



1-1 



2a A\ 
■ In 



3tt 



m/ 



327r 2 9m 2 2 



1024tt m 4 a. 



2a s // B 
1 In 

37T 3m 



1-1- 



37r 3m / 



(40) 



(41) 
(42) 

(43) 
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Adding both contributions, we obtain 

_ ^y(vac) _|_ £,(mat) 



3e 9 y2 

32^ 



(44) 



c 2 = cl vac; + c 2 

,4 



e 5 



\-Z$ (jT 3 - fi m 3 ) , (45) 



1024vr m 4 a_ 
where 

fi = l_^£l n A( =Z2 ) (7 m = i_^£i n ^£. ( 46 ) 

37r m 37r 3m 

4 Effective potential 

In this section, we study the structure of the effective potential V((jl)). The "gener- 
ating function" W(B 2 ) is related to (/I'), Eq. (27), through 

w(^)=/ o V>^=^ 2 + ^(^ 2 ) 2 (47) 

with Ci and C 2 as in Eqs. (44) and (45), respectively. V ({£?)) is defined by the 
following Legendre transformation on W (-B 2 ): 

K«/r» = BB-lt'(B') , (48) 

= B, (49) 

where B {= B ((/?'))) ^ s related to (/T) through Eq. (27). 

When one is interested in V ((/*')) with small (//'), it is sufficient to obtain B ((//)) 
in powers of (/7/). For obtaining more accurate information, however, we solve Eq. 
(27) exactly. 

As will be seen below, in the region of our interest, C 2 > 0. For convenience, let 
us introduce the dimensionless quantities, 



3 / M 
M = 2i\C^^' ) = M \M\ {M -° h (50) 



V(M) ee (<//)). ( 5 2) 
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Eqs. (27), (48) and (49) turn out, in respective order, to 
B 3 (M) + 3e(d)B(M) - 2M = , 
V(M) = 2 -MB(M) - ^e(d)B 2 (M) - ^B\M) , 

dM 9 v ; 
where e(d) = Ci/|Ci|. 

Solving Eq. (53), we obtain the following expressions for B(M). 

d > 

B(M) = (Vl + M 2 + M) 1/3 - (Vl + M 2 - M) 1/3 . 

From Eqs. (55) and (56), we see that 

dV{M) 



dM 
dV(M) 

dM 



> for M > , 
= at M = . 



d < 



(1) M > 1 

Solution (si): 

jB(m) = (m + Vm 2 - l) 1/3 + (m - Vm 2 - l) 1/3 

dV(M) 



dM 

(2) M < 1 

Solution (s2): 



> 0. 



Solution (s3): 



*w„,n ^ /arccosM\ 
tf(M) = 2cos(^ J : 

g) 

47r + arccos M 



B(M) = 2cos( 
dV(M) 



dM 
dV(M) 



3 

< for < M < 1 , 



= at M = . 
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Solution (s4): 

'2n + arccosM 



„/,,>. n /27T + arccosM\ 
B(M) = 2cos(^ J 

dV{M) 



AM < °- 

In Fig. 2, we display V(M) in the case of C x > 0. The curve "0" is V(M) with 
Eq. (56) for B(M). We see that, when C\ > 0, the quark matter is in the normal 
phase. 

In Fig. 3(a), we depict V(M) for C\ < 0. The curves (si), (s2), (s3) and (s4) show 
V(M) corresponding, in respective order, to the solutions (si), (s2), (s3) and (s4) for 
B(M). The actual effective potential V(M) is as in Fig. 3(b). From this figure, we 
see that, when C l < 0, V(M) is minimum at M = 1 or (p') = y / 4|Ci| 3 /(27C 2 ), which 
means that the spontaneous magnetization occurs. 

Here we discuss the effect of the higher-order (in B) contribution to (p') in Eq. 
(27), 

(pf) = [d + C 2 B 2 + C 3 (B 2 ) 2 } B. (57) 
Eqs. (53) and (54) turn out to be 

B 3 (M) + 3e(d)B(M) + ^j^B 5 (M) - 2M = , (58) 

V(M) = 2 -MB(M) - \e{C x )&{M) - ±-B\M) - ^f^(M) . (59) 



From these equations, we see first of all that, near the phase transition point, C\ ~ 0, 
the effects of the "C3 term" is very small. We will compute V(M) for the following 
two choices for \Ci\Cz/C 2 : 

= -2, +2. 



Cl 

The curve "+" ("-") in Fig. 2 depicts V(M) with \d\CsfC 2 = +2 (-2). We see that 
the effect of the "C 3 term" in Eq. (57) on V(M) is small and the system is in the 
normal phase. Fig. 3(c) depicts V(M) in the case of C\ < 0. For |Ci|C 3 /Cf = +2, the 
behavior of V(M) is qualitatively the same as that for C3 = and the spontaneous 



magnetization takes place with (p 1 ) ~ 0.84y4|Ci| 3 /(27C2). On the other hand, for 
\C\\C 3 jC\ = —2, V(M) decreases monotonically as M increases. If we take this result 
at face value, the system is in the ferromagnetic phase. 
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Recalling the fact that C\ is a function of A and fi B , we summarize the above 
observations: 

• In the region of /i B where C*i(A, n B ) > 0, the system is in the normal phase. 

• In the region where Ci(A,fi B ) < 0, the system is in the ferromagnetic phase. 

• At 

/Ab = /Abo where C \(A, fi Bc ) = 0, a magnetic phase transition between the 
above two phases tales place. 

5 Numerical Analysis 

In this section, on the basis of the observation in Sec. 4, we study the magnetic 
property of the quark matter through the numerical analysis in the region 0.6 GeV 
< < 4 GeV. Let start with two observations. 

• As the renormalization scale hr (for the running coupling constant a <,(//#) and 
the running mass m(/XR)), we choose fiR = fi B . There is an ambiguity on the 
choice of hr. For computing the pressure of the quark-gluon plasma (QGP) 
at temperature T and fi B = 0, /ir ~ 2tcT is chosen in the literature (see, e.g., 
Ref. 9)). The so-called thermal mass of a zero mass quark in QGP with finite 
T and fi B is m t h oc T 2 + n B /(37T 2 ). These observations suggest the choice, 
li R = 2vr[T 2 + /il/^Tr 2 )] 1 / 2 ^4 2tt x (fji B /(y/3ir)) ~ 1.15/i B , which is not very 
different from /i B . 

• Following Ref. 10), we identify the effective potential V computed above with 
cutoff A with the one renormalized at the renormalization scale jiR = A in QCD. 

On the basis of these observations, as a crude approximation, we choose A = /jl b . 
Then, from Eqs. (45) and (46), C2{^b) > 0, provided that Q((J> B ) > 0, which is the 
case in the region of our interest. 

For the QCD a s (= g 2 /(47r)), we take the running one a s ([i B ) in the MS scheme 
with a s (fi B = 1.2 GeV) = 0.39. In Fig. 4, a s (fi B ) is depicted in the relatively small 
fi B region, 600 MeV < fi B < 700 MeV. For the quark masses, we take 1-loop running 
masses 
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with 

m u (p R = 2 GeV) = 3 MeV/c 2 , m d {p R = 2 GeV) = 6.8 MeV/c 2 , 
m s (p R = 2 GeV) = 118 MeV/c 2 , m c (p R = 1.2 GeV) = 1.2 GeV/c 2 . 

rif is the number of quark flavors and the quark with mass m q is counted for m q < p R 
(— A*b)- F° r the QED fine structure constant a, we take the 1-loop running one. For 
computing it, we have used the above-mentioned quark masses and, for lepton masses, 
we have used m e = 0.5 MeV/c 2 , m M = 106 MeV/c 2 and m T = 1.78 GeV/c 2 . 

The relation between the baryon-number density p and the baryonic chemical 
potential p R is 

i / \ / 2 \ 3/2 

p{p B ) = 3 x 2 x — X; (9 ^ - m.fe) J ^ - m 2 (/i B )j , (60) 

where "3" is the number of colors and "2" comes from the spin degree of freedom. 
We do not take into account the QCD correction to the relation (60) except that 
the running masses have been used. The nuclear density at the center of a nuclei is 
Po ~ 0.17 fm -3 . The density of neutron stars are in the range, p = 10~ 3 po ~ 10po- 
For a guide of eyes, the relation between p/p and p B is shown in Figs. 5(a) and (b), 
and the ratio pi/p (i = u, d, s, c) is shown in Fig. 5(c). 

Quarks 

As mentioned above, C^/Ab) > in the region of interest. Numerical analysis shows 
that 

{Ci(fi B ) < for p B < psc , 
Ci(p B ) = for p B = Pbc , 
Ci(hb) > for pb > Hbc, 
with psc — 641, 625 and 693 MeV for, in respective order, u, d and s quarks. In Fig. 
6, Ci(ps) is depicted against ps for 600 MeV < ps < 700 MeV, where the coefficient 
d\ = 1, 50 and 2 x 10 4 for u, d and s quarks, respectively. Vertical dashed lines 
indicate the phase transition points. We see that, at high-density region, p B > I^Bc, 
C\ > and the quarks are in the normal phase. As the density is lowered, the 
quarks undergo a magnetic phase transition (Ci = 0), at the critical density, into 
a ferromagnetic phase {C\ < 0). We see from Fig. 3(b) (also Fig. 3(c)) that, in 
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the ferromagnetic phase, there arises spontaneous magnetization whose magnitude is 



K/0| = V /4|C 1 |-V(27C 2 ) M ~ X /4|C 1 |V(27C 2 ) (= (p'>) (cf. Eq. (50)). We depict (pi') 
in Fig. 7, where <i 2 = 10" 4 , 2 x and 0.1 for u, d, and s quarks respectively. 

Ci is the quantity that determines the phase of the system. Eq. (40) with A = p B 
tells us that the vacuum-sector contribution to C\ is negative, c[ va ' c ^ < 0, in the 
region of our interest. While, from Eq. (42), C[ ma ' t \p B ) is positive and an increasing 
function of p B . The effect of the QCD interaction does not change the signs of C{ vac) 
and C{ mat) . At low densities Cf nat) is so small that d = Cf ac) + C{ mat) < 0. As p B 
is increased (7j mat - ) increases and we arrive at C\ = at p B = p Bc . As p B further 
increases, C\ turns out to be positive. To understand the behavior of C[ mat \p B ), let 
us recall that the hamiltonian density describing the interaction of the system with 
the external magnetic field is 7i m t = —p! ■ B. This means that the configuration 
p! || B is energetically favorable, which reflects on c{; mat ) _ Q > 0, Eq. (30) or (42) 
(cf. Eq. (27)). 

Several observations are in order here. 

• In spite of the fact that m u < < m s , we have obtained the curious result, 
{fJ-Bc)d < {^Bc)u < {f^Bc)s- This means that, as fi B is lowered from the normal- 
phase region, the heaviest s quarks are magnetized first, then follows the lightest 
u quarks and finally d quarks start are magnetized. 

• The magnitude of the magnetization is larger for the lighter quark. This is a 
reflection of the fact that jl' oc 1/m. 

• The phase transition points [i Bc ~ 641,625 and 693 MeV for, in respective 
order, u, d and s quarks correspond to the densities p ~ 1.62p , l-47po and 
2.20p with p the nuclear density. These values for p are near to the nuclear 
density p an d then our results may be relevant to the quark star. 

• As seen from Fig. 4, a s (p B ) is not very small at 600 MeV < p B < 700 MeV. 
Then, before arriving at definite conclusion, higher-order (in a s ) contributions 
should be taken into account. 2 



2 In this relation, we recall the fact that, for a dense electron system, the higher-order contributions 
substantially changes [11] the magnetic phase diagram when compared to the calculations using the 
random-phase approximation. 
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Here we mention the case with the order parameter p = 2p' (see Appendix A). It 
can readily be seen that C x = 4Ci(Eq. (44))/4 and C 2 = lQC 2 {Eq. (45)). The phase 
transition point p Bc does not change. 

Leptons 

Our primary interest is to study the magnetic property of the quark matter, like quark 
star. Such a system is electrically neutral, so that there exists leptons. From the 
charge neutrality, we obtain the relation between the baryonic chemical potential p B 
and the leptonic chemical potential pl being conjugate to the lepton number. We see 
that, for 0.6 GeV < p B < 4 GeV, p L covers the range 174 MeV < p L < 969 MeV. In 
Fig. 8(a), we depict the relation between them in the region 0.6 GeV < p B < 3 GeV, 
and, in Fig. 8(b), the ratio pi/p (i = e,p) is shown in the region 174 MeV < p L < 
300 MeV. An abrupt increase in p L at p B ^ 2.9 GeV is due to the opening of the c 
quark channel, m c (p B ^ 2.9 GeV) ~ 2.9/3 GeV/c 2 . 

For the QED coupling constant, we take a = 1/136 and, for the lepton masses, 
we take the values cited at the beginning of Sec. 5. Numerical analysis shows that, 
in the region 174 MeV < pl < 969 MeV, the leptons are in the normal phase. 

6 Summary and Outlook 

In this paper, we have studied the magnetic property of the quark matter through 
evaluating the effective potential V( (p 1 )) for the magnetic moment p! of a quark. We 
have found that, at low densities p < p c , the quarks are in the ferromagnetic phase. 
For the u, d and s quarks, p c ~ 1.62p , l-47po and 2.20po (po the nuclear density), 
respectively. The phase transition occurs at p = p c and, at high densities p > p c , the 
quarks turn out to be in the normal phase. 

We have also studied the case of leptons in the region 174 MeV < pl < 969 MeV 
(Pl the leptonic chemical potential), which corresponds to 0.6 GeV < p B < 4 GeV. 
We have found that both e and p are in the normal phase. 

We list the several points which are left for future study. 

1. When more than two kinds of particles are in ferromagnetic phases, for deter- 
mining the relative direction(s) of their magnetizations, interactions between 
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them should be taken into account. 

2. Study of the quark matter at finite temperature and density. 

3. Analysis in the region 'CSC in Fig. 1, where several phases coexist. 

4. Undoing the assumption made after Eq. (1) in Appendix A. 

5. Improvements of the approximation. 

a) Undoing the leading logarithmic approximation (1). 

b) Resummations for the (quark and gluon) propagators and the quark-gluon 

vertex. 

c) Computation of the higher-order contributions. This is particularly impor- 

tant for studying the ferromagnetic phase, since cx s ([/,b) is not very small 
there. 

d) Resummation of, e.g., the ladder diagrams. From thus obtained effective 

potential and the renormalization-group formula, one can expect to obtain 
(approximately) A-independent result. (For the case of chiral and di quark 
condensations, see, e.g., Ref. 10).) 

6. Computation of the effective action, which allows us to study the dynamical 
aspect of the system. 

7. Study of nonequilibrium quark matter, which allows one to deal with the space- 
time evolution of the system under a given initial data. 
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Appendix A Order parameter 



We start with the action that describes the interaction between a quark (with charge 
e q ) and uniform external magnetic field B: 



I 



d xL 



cxt 



d 4 xip^f ■ Aip 



= y/ d 4 x^-{B xr>, 



where tp and ip are the quark fields in the interaction representation. Here the color 
index is suppressed. Using the Dirac equation, we have 

= ip(x) I —% — m J 7 • A(x)ip(x) + 4>{x)~f ■ A(x) (i0 — mj ip(x) , 



which yields 



2771-07 ' Aty = ip 



-iA-V + A%a iik + i di A- <9 M cr^A 



1> 



— * — * 

Using A = B x f/2, we finally obtain 



J d 4 xC, 



ext 



= [d 4 xip \-2iB • (f x V) + a ki e ijk B j ] 

Am J L J 

<„ I ., . - ^ 



= ^- I d 4 x 

2m 



-ifxV + a 



We take a statistical ensemble of the systems, each of which has vanishing to- 
tal angular momentum. We further assume that, in each system, the total angular 
momentum of quarks and that of antiquarks vanish separately. Incidentally, the an- 
tiquarks exist as the virtual particles, through repetition of the gg-pair productions 
(from the vacuum) and the pair annihilations. Then, Eq. (1) turns out to 



/ 



d xL 



cxt 



/' 



d A x ■ B , 



As the order parameter, we adopt jl' . 

It is to be noted that if we ignore the contribution from the orbital angular mo- 
mentum in Eq. (1), we obtain Jl (= 2/2') for the order parameter. 
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Fig. 1. Sketch of the QCD phase diagram. 'T' is the temperature and £ /xb' is the 
baryonic chemical potential. 'Hadron' stands for the hadronic phase, 'QGP' for 
the quark-gluon plasma phase and 'CSC for several phases that include the 
color super conducting phases. 

Fig. 2. V(M) in the case of C x > against M. The curve with "0 " is the V(M) 
with Eq. (56) for B(M). The curves with "+ " and "- " is V(M) in Eq. (59) 
with |C1|C 3 /C| = +2 and —2, respectively. 

Fig. 3. V(M) in the case of C\ < 0. (a) The curves (si), (s2), (s3) and (s4) are 
V(M) with, in respective order, the solutions (si), (s2), (s3) and (s4) for B(M). 
(b) Actual effective potential, (c) The curve with "0 " is as in Fig. (b) and the 
curves with "+ " and "- " are V(M) in Eq. (59) with |Cl|C 3 /Cf = +2 and -2, 
respectively. 

Fig. 4. a s {fi B ) against 600 MeV < p B < 700 MeV. 

Fig. 5. (a) The ratio p/ p against 0.6 GeV < p B < 1-6 GeV. (b) Same as (a) for 
0.6GeV < p B < 4 GeV. (c) The ratio pij p (i = u, d, s, c) against p B - 

Fig. 6. Ci(// B ) against 600 MeV < p B < 700 MeV. d x = 1,50 and 2 x 10 4 for u, d 
and s quarks. 

Fig. 7. \{p') \ = v /4|Ci| 3 /(27C 2 ) against 600 MeV < p B < 700 MeV. d 2 = 10" 4 , 2 x 
10~ 3 and 0.1 for u, d and s quarks. 

Fig. 8. (a) Relation between p L and p B . (b) The ratio pi/p {i = e,p) against p^. 
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